In this paper, we are going to approximate the Dynamical system of the SIR mathematical model, numerically. The expression for basic reproduction number R0 is obtained which plays main role in the local and global stability. On the other hand, if R0 > 1 then the equilibrium point is unstable. Our numerical simulations show that our proposed NSFD method provides reliable and consistent, positive and converging results at all step sizes as compared to the traditional RK-4 method and Euler method which diverge at large step sizes and produce negative and unstable results with large oscillations. The local stability of the equilibrium point is proved by Jacobean method,however, Schur-Cohn conditions are used to discuss the local stability for the discrete NSFD scheme. To prove the global stability of equilibria, the Enatsu criterion and Lyapunov function are used. This paper presents theoretical findings and numerical simulations that can be used as a powerful instrument for forecasting the spread of added infectious illnesses.